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D Abstract 

In this paper we study the properties of the phase diagram of a simple extra dimensional 

,,^ , model on the lattice at finite temperature. We consider the five-dimensional pure gauge 

^ '. 
cn ' abelian model with anisotropic couplings which at zero temperature exhibits a new inter- 

\o . 

^ , esting phase, the layer phase. This phase is characterized by a massless photon living on the 

four dimensional subspace and confinement along the extra dimension. We show that, as 



long as the temperature takes a non zero value the aforementioned layer phase disappears. 



OO , It would be equivalent to assume that at finite temperature the higher-dimensional lattice 

model loses any feature of the layered structure due to the deconfinement which opens up 



the interactions between the three-dimensional subspaces at finite temperature. 
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1 Introduction 

The original idea of Fu and Nielsen in the mid-eighties of the last century consisted in proposing 
a new way of dimensional reduction through higher dimensional lattice models with anisotropic 
couplings [HI. Since then in a series of papers the phase diagram of the higher-dimensional lat- 
tice models was studied by mean field and Monte carlo methods. Besides that a mechanism 
of producing the anisotropic couplings was proposed invoking a Randall- Sundrum space-time 
metric to the (continuum) higher dimensional models [l2l[3l. The lattice model with anisotropic 
couplings which came up could help in understanding the localisation of the gauge interaction 
on the brane; the idea was that the usual four-dimensional space-time is embedded in a higher 
dimensional bulk in which the extra dimensions are subject to confinement. Indeed the assumed 
strong coupling dynamics along the extra dimension requires a non-perturbative study. The nu- 
merical study on the lattice has verified the prediction of a new phase (layer phase) proposed by 
Fu and Nielsen who used mean field methods. In this new phase we established the existence of 
a massless photon on the four-dimensional subspace while at the same time the extra dimension 
is confined [4J . This confinement is responsible for the fact that the higher-dimensional space 
is layered-like and the interactions are confined in the four-dimensional space-time slices. 

It is worthwhile to mention that in the Dvali-Shifman model a similar way of thinking has 
been used in order to achieve a localization mechanism on a brane: the assumption of con- 
finement (of a non-abelian nature) along one of the dimensions which limits the dynamics of 
the model in a subspace with one dimension less 0. It is possible to generate gravitationally 
this localization mechanism for the trapping of the charged fields, under a non-abelian gauge 
field, on a 3d submanifold (brane) using the non-minimal coupling of gravity with a scalar field 
O. The result is a spontaneously broken phase on the brane (higgs phase) and a confining 
(symmetric phase) in the transverse directions (bulk). 

In this work we consider the following exercise: we assume a five-dimensional U(l) lattice 
model with anisotropic couplings at finite temperature. Our intention is to discover the fate of 
the layer phase as we switch on the temperature to non zero values. By means of numerical 
methods we study the phase diagram and our main result is that for T 7^ the layer phase 
becomes a deconfined phase with new properties. In other words, the confining extra dimension 
that is detected at T = 0, becomes deconfined for non zero values of the temperature and the 
system is lacking the four-dimensional layered structure, [j 

The previous result is suggestive for the behavior of a multidimensional anisotropic gauge- 
higgs model at finite temperature a It is well known that at high temperature the higgs phase 
turns into a symmetric phase. We expect then that the layer higgs phase disappears at this 



'Our results contradict the prediction based on a Variational Cumulant Expansion by the authors of ref. |7 1 for 
which the layer phase of the anisotropic lattice gauge models persists at high temperature. 
^For the analysis of anisotropic gauge-higgs models at zero temperature see the refs. |l8]|9l. 



temperature. The new phase is probably a muhidimensional high temperature symmetric phase: 
a hot multidimensional world in the "quark-gluon plasma" phase. In this rather hypothetical 
scenario the Universe starts as a hot multidimensional system that cools down, passes through 
a series of phase transitions and ends up to a brane Universe at zero temperature. 

To explore the phase diagram of the anisotropic 5D U(l) gauge model we have to under- 
stand first the phase structure of the 4D U(l) gauge model at high temperature. The study of the 
phase structure of lattice electrodynamics in three and four dimensions at zero temperature, us- 
ing the topological excitations of the theory (monopoles) was first performed by the authors of 
ref . [iOj . Computer simulations for the 3D compact QED at finite temperature were performed 
in ref. ifTTI for the deconfinement transition from the monopole anti-monopole point of view. In 
three dimensions and zero temperature the theory is confining for all values of the coupling con- 
stant and a monopole and anti-monopole plasma is responsible for the permanent confinement 
of oppositely charged particles. For non-zero temperature the binding of monopoles and the for- 
mation of magnetic dipoles lead to loss of confinement. The dipole plasma can not sufficiently 
screen the field created by electric currents and the screening mass vanishes. 

In four dimensions the point-like topological excitations become one-dimensional objects 
(strings of monopole currents). Again, in the zero-temperature case and for small values of 
the coupling constant P there is a large number of monopole loops winding around the system 
and causing disorder. As a consequence, if an external field is applied it will be shielded after 
a small penetration. For large values of P, on the other hand, the situation is different. A 
long distance penetration of the external field is observed, accompanied by the renormalization 
of the magnetic charges due to the monopole currents. The case of 4D compact U(l) gauge 
theory at finite temperature was studied separately in fT2l| and fTSl . The authors of reference 
lfT2l reported a second order phase transition to a Coulomb phase for Lt > 4, with critical 
exponents consistent with 3d Gaussian values and no obvious dependence on L^. A different 
picture emerged in ref. |fT3l, where among other things the disappearance of the Coulomb phase 
for all values of T{= x") 7^ was predicted. Instead of a Coulomb phase we are left with a 
spatial confining - temporal Coulomb phase for all temperatures. 

In this paper, we are not going to present a detailed study of the nature of the phase tran- 
sitions; however, some of our findings seem to indicate the absence of a Coulomb phase for 
all temperatures different from zero for the case of four dimensional compact QED (as long 
as the condition Lt « Lg is satisfied). We go one step further and examine a finite temper- 
ature scenario in five dimensions through the anisotropic U(l) gauge model with couplings (3 
and (3 . The connection of this model with the brane model scenarios makes it an ideal can- 
didate for the study of the brane models in the non-zero temperature case. We are mostly 
interested to discover if some of the most promising characteristics of this model survive in the 
high-temperature regime. In what follows we will try to give a brief summary of our findings 



through the description of the limiting cases of our model. 

For /?' = we obtain the four dimensional QED at finite temperature. From the study of 
the system with volume V^d = Lt x L^ and L^ = 2, 4 and 6 we come to the conclusion that, 
although phase transitions seem to appear for finite L^, they are actually finite-size effects and 
disappear in the limit Lg ^ oo. For all values of L^ 7^ the Coulomb phase gives its place to 
a temporal Coulomb - spatial confining phase (deconfining phase) lfT4ll . The Coulomb phase is 
recovered only at T = 0. 

For P ^ and Lj = 2 we have a five dimensional, anisotropic model in a high temperature 
state. The zero temperature model has been already studied and it is characterized by three 
distinct phases [|3l|4l[T5l[T6l. A five dimensional confining phase, a 5D Coulomb phase and the 
layer phase where the system is confining along the fifth direction while, along the four remain- 
ing directions, it exhibits the Coulomb behaviour with a massless photon. These characteristics 
change when the temperature becomes non-zero. We observed the replacement of the layer 
phase by a deconfining phase, due to the same mechanism responsible for the disappearance 
of the Coulomb phase in four dimensions. The behavior of the system in the time direction is 
coulombic and confining in the remaining four directions. 

Our work is organized as follows. In section 2 we present the action of the model and the 
observables, the helicity modulus and the Polyakov line, that we use in order to characterize the 
phase diagram of the model. In section 3 we analyze the system in the three limiting cases: 5D 
anisotropic for T = 0, /3 =0 and Lt = 1. Finally in section 4 we present thej)hase diagram 
for the 5D anisotropic model at finite temperature and in particular for Lt = 2 
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2 The model 

2.1 Definition 

The five dimensional anisotropic U(l) gauge model with two couplings, [3 and , at finite 
temperature is defined as : 

x,l<^<i/<3 3:,l</i<3 

/3' Y. {^~Re{U,,{x))+[3'Y,(l~Re{Ut,{x)) (1) 

a;,l</i<3 X 



^We recall that T = l/Lt in lattice units. 



where 

Uf,u{x) = U^{x)U^{x + asfi)Ul{x + asi>)Ul{x) 

U^t (x) = U^, (x) Ut {x + as ft) Ul {x + ati) U} {x) 
Uf^bix) = U^{x)U5{x + asfi)Ul{x + a55)Ul{x) 
Ut5{x) = Ut{x)U5{x + ati)U^t{x + a55)Ul{x) 

are the plaquette variables defined on the 4d-subspaces {(//, u = 1,2,3)- t} and planes con- 
taining an extra, fifth dimension (X5). With an obvious noatation we call these plaquettes as 
Ps,PsuP's,^ndP;,. 
The link variables are defined as 

U^{x) = e'^-("), Ut{x) = e'''*^") and U^ix) = e'^^^") 

Let us make the notation clear. The action is defined in an Euclidean lattice volume, namely 
V = Lt X L^ X L5 in lattice units. Lt is the compactified temporal dimension which is related 
to the temperature through the relationship 

T = -^ (2) 

Uat 

We denote with at the lattice spacing, Lt is an integer number, -Ls=i 2,3 are the usual "infinite" 
space dimensions and finally L5 is an extra, fifth dimension, which we consider to be infinite and 
equal to Lg. We assume periodic bourdary conditions for the U(l) gauge field in all directions. 
The proclaimed anisotropy of the model has nothing to do with the "time" direction. In this 
model the lattice spacings as,at are equal. The anisotropy is introduced through the interaction 
along the extra direction. So, we have ag = at = a and a^ ^ a where a^ is the lattice spacing 
related to the extra dimension. 

In our model the gauge couplings /3 and P' are generally independent from each other and 
the coordinates. The lattice spacing is determined from the value of the couplings P and P . In 
some cases we can have a coordinate dependence and it is possible to relate them with extra 
fields, as in the brane model [l2l|3l[T7i|. In terms of the continuum fields the link angles, 9m, can 
be written as : 

Om{x) = aMAM{x) 

where Am{x) are the gauge potentials [[31 [8l and with M we denote M = (t, /i, 5). In the naive 
continuum limit (a, a^ -^ 0) we define: 

/3 = ^ and /3' 



9I gla^ 

where g^ is the bare five-dimensional coupling constant for the gauge field. The resulting con- 
tinuum action takes the standard form: 

J 95 



Note that g1 has dimensions of length and is related to a characteristic scale for five dimensions. 
The previous expression does not exhibit any anisotropy at all. However, the results that we 
present below indicate that the anisotropy may survive in the continuum limit. 

2.2 Observables 

We now proceed to the introduction of the observables, i.e. the gauge invariant quantities which 
are used for the study of the model. 

2.2.1 The helicity modulus 

Among the quantities used to distinguish the various phases and the respective phase transitions 
in a statistical model the ones that attract the most attention are the so called order parameters. 
Their great significance comes from the fact that they display completely different behavior 
between the various phases. Their "thermal average" is zero on the one side of the transition 
and moves away from zero on the other side. For the case of a confining-Coulomb transition 
a quantity with the properties of an order parameter is the helicity modulus (h.m). It was first 
introduced in the context of lattice gauge theories by P.de Forcrand and M. Vettorazzo and it 
characterizes the responce of a system to an external flux. It is zero in a confining phase and 
nonzero in a coulombic one[13]. 

Let us consider our five dimensional system with (L^, L^, Lp, L^, L5) and let us choose a partic- 
ular orientation, for example, {fi, v). Through the remaining orthogonal directions it is defined 
a stack of Lp x Lt x L5 plaquettes parallel to the (ji, u) orientation. In order to study the re- 
sponse of the system to an external static field we assume the presence of an external flux $ 
through this stack of plaquettes. By a suitable choice of variable transformations we can spread 
the flux homogeneously over the parallel planes. In other words, we can add the constant value 
of $p = -j^j^ to each of the plaquettes of the given (/x, z/) orientation. Also we can impose 
an external flux by changing the boundary links using twisted boundary conditions instead of 
using pure periodic ifTBllTSl . The partition function, in the presence of the external flux, is: 

Z($)= [ DOe-^^'''''^ (3) 



{^u)planes ^ {fiu)planes 

-P ^ cos{9^tix)) - P' ^ cos(6'^5(x)) 

-/3'J2^os{et,{x)) (4) 



where Yj(nu)pianes ^^ the sum over the plaquettes of the given orientation (/i, u), containing the 
flux and Yj(ju7)pianes ^^s Complement, consisting of all the plaquettes that remained unchanged 
(plaquettes belonging to the other planes). 
^From the partition function we can obtain the flux dependent free energy 

F($) = - ln(Z($)) = - In r /" DOe^^^^'^^A (5) 

An important observation is that the partition function Z{^) of equation (3) and hence the flux 
free energy is clearly 2n periodic . So, the extra flux we impose on the system is defined only 
mod(27r). 

In the confining phase the flux free energy -F($) is constant in the thermodynamic limit because 
the correlation length and the effect of the external flux through the twisted boundary links is 
exponentially decreasing. On the contrary, in the Coulomb phase we have an infinite correlation 
length, so the influence of the twisted boundary conditions is extended to the full extent of the 
system. As a result we have a nontrivial dependence of -F($) by the external flux $. 
The helicity modulus is defined as 






5$2 



(6) 

$=0 



and it gives a measure of the curvature of the free energy profile around $ = 0. From the above 
equation and for various choices with respect to the orientation, due to the anisotropy of the 
model, we have: 



hs{(3) = jj^, ((T{Pcos{e,^))\ - /(^(/3sin(V)))A j (7) 



{L^Lt 




hm = — ^ f ( j2(^cos{e,,)) ) - ( Q2{f3M0,t))r ) ] (8) 



hs,{P' 



(L 



i-p ( r£0 cos(^,5))\ - l{J2(^' MO,,))A ] (9) 



ht,i0) = jT^r^ I ( Ti^' cos(^,5)) ) - ( :i0 MOt,))r ) I (10) 



{LtL 



5) 



P' I \ P' 

to to 



Now, consider for the moment the classical limit (P -^ oo) for the action (4) where all the 
fluctuations are suppressed. In this limit the flux is distributed equally over all the plaquettes 
of each plane and it does not change as we cross the parallel planes. If we expand the classical 



action in powers of the flux, since in the thermodynamic limit the quantity j^^ is always small 
, we find: 

1 n^2 ^5D , ^_^ ^ j^ /^N 77 rn\ ^ Qif.2 ^5D 



5'c,a.,cai($) = -/3^ /r T N2 + ^onstant =^ F,„,,„„($) - F,„,,„„(0) = -/3$ 



2r^ ■*- / J J \o *^»jiiuvt.*.iiv / J. classical \ -*- / -^ classical \'"'/ o ' / 7" 7" \ 9 

where Vs^ = L^^L^LpLtL^ is the 5D lattice volume. 

The above expression for the free energy, F, holds all the way up to the phase transition, 
where fluctuations are present, if one only replaces the bare coupling by a renormalized cou- 
pling, (3 -^ (3r{I3) (for details see lfT8l[T3l): 

^From the Eqs. (6) and (1 1) we have for the "spatial" h.m 

K{(3)^(3rU (12) 

and following the same steps, we can get the scaling relations for the remaining quantities: 

htW) ~ /?ij/ (13) 

h,,{(3') ~ (3'j,Lt (14) 

ht5iP')r-PR^ (15) 

Although the arguments presented above are based mainly on the classical approach, this is 
indeed the case in the Coulomb phase and the helicity moduli applied for the five dimensional 
system behave exactly as the above equations predict. 

2.2.2 Polyakov loop (or Wilson line) 

For the evaluation of the potential between a static quark-antiquark pair at zero temperature, 
the study of the ground state expectation value of the Wilson loop for large Euclidean times is 
needed. When the temperature is non zero (Lt << Lg as opposed to the former case) the same 
information is obtained by using a different quantity which is the Polyakov loop or the Wilson 
line. It consists of the product of link variables along topologically non-trivial loops, winding 
around the time direction due to periodic conditions. 

Lt 
nt=l 

^* = ZJ E ^*(^'^5) (16) 

* (n,n5) 



where {n}eZ^ denotes a lattice site. 

Physically, the expectation value of the Polyakov loop determines the free energy of a system 

with a single, static heavy quark, measured relative to the vacuum: 



(|P,|) = e~^^'^ (17) 

where |-P| is the absolute value of P and < ■ ■ ■ > the statistical average value evaluated using 
the action of equation (1). The above relation holds even in the presence of finite mass quarks 
coupled to the gauge potential with the only difference that in that case the expectation value 
has to be calculated using the finite temperature action that includes the dynamical fermions. 

The Polyakov loop is somewhat the world line of a static quark in a Wilson loop and that 
suggests that the free energy of a quark-antiquark pair located at (ni, nsj and (tt^, ^53) respec- 
tively is given by the correlation function of two such loops, with bases at the aforementioned 
points and having opposite orientations. Consequently we have: 

(P,,(ni,n5j4(n2,n5j) = e-i*^,,-(Ki.«5j;{nl,n5,}) (ig) 

keeping Lt constant. 

For large distance separation of the quark-antiquark pair and assuming that the correlation func- 
tions satisfy clustering, the above expression reduces to 

Pt,{ni,n5jPl{n2,n5^)\ -^ \ {Pt) |^ for \R\ -^ 00 (where P = {nl.nsj - {nsj^sj) 

(19) 
which is just the self-energy of two isolated quarks. 

In the confinement phase the correlation function of the Polyakov loop decays exponentially for 
large distances: 

'Pt{0)PhR)) ~ e-^*'^l^l (20) 



giving a linear potential with string tension a and Fgg ~ a\R\. 

The flux free energy Fgg increases, in general, for large separation of the quarks in the confining 
phase, giving eventually < \Pt\ >= and Fgg = 00 in the thermodynamic limit. We interpret 
< \Pt\ >= as a signal for confinement. If we have < \Pt\ >^ 0, then the free energy of the 
static quark-antiquark pair tends to a constant value, for large separation of the heavy charges, 
as shown in Eqs. (18) and (19), and this is a signal for deconfinement. In other words, the ex- 
pectation value of the temporal Polyakov loop serves as an order parameter in finite temperature 
gauge theories. 



3 Three limitting cases 
3.1 The zero temperature case 

The five dimensional anisotropic U(l) gauge model, at zero temperature, was first introduced by 
Fu and Nielsen [[H as an attempt to offer an alternative way to achieve dimensional reduction. 
Since then many numerical investigations of the model have been made lfT5l[T6ll . As we have 
already noted in the Introduction, the interest in this comes from the fact that the anisotropy of 
the model produces a new phase, the so called layer phase, which can serve as a mechanism 
for gauge field localization on a brane. We can induce this anisotropy to the gauge coupling 
using, for example, the Randall- Sundrum metric background in five dimensions. The effect of 
the warp factor from the RS background or a general anti-de- Sitter (AdS^) background on the 
U(l) gauge theory is to provide the gauge theory with a different gauge coupling in the fifth 
directiondd). 

In Fig. 1 we present the phase diagram of the theory. It consists of three distinct phases. For 
large values of /? and /3 the model lies in a Coulomb phase (C) on the 5-D space. Now, if one 
keeps /3 constant and bigger than one and at the same time lets decrease /?', one will eventually 
come across the new phase, the layer phase (L), where the forces in four dimensions will still be 
Coulomb-like but in the fifth dimension the confinement is present. For small values of P and /?' 
the force is confining in all five directions and the corresponding phase is the Strong phase (S). 
The properties of the three phases can become more transparent in terms of two test charges. 



The phase diagram for zero temperature 



CO. 




Figure 1: The phase diagram for the 5D anisotropic U(l) gauge model at zero temperature. 
Three phases are present: Strong confining phase, 5D Coulomb phase and the Layer phase. 
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In the Coulomb phase the force between two heavy charges is 5D Coulomb-like, and becomes 
the exact five dimensional Coulomb law in the diagonal line, defined hy (3 = for which 
no anisotropy appears (for a numerical investigation see flU). The completely opposite picture 
emerges in the Strong phase. There the force is confining in all five directions giving infinite 
energy for the separation of the test charges in any direction. Now, two test charges in the 
layer phase will experience a Coulomb force in the four dimensional layers, with the coupling 
given by the four-dimensional coupling (5; there are strong indications of the similarity with the 
usual 4D Coulomb law (see BU for details), while along the fifth direction the test charges will 
experience a strong force as the corresponding coupling /5 takes small values. This means that 
charged particles in the layer phase will mainly run only along a layer since in an attempt to 
leave the layer in which they belong they will be driven back by a linear potential, analogous to 
the one responsible for the quark confinement. This is the mechanism on which the gauge field 
localization scheme is based. 

Now we would like to sketch the three phases in terms of the helicity modulus. In the zero 
temperature case (L^ = Lg = L5) we are left with only two possible choices. Instead of the 
Eqs. (7)-(10) we have: 

hsm = j^^, ({TWcos{e,^))\ - /($^(/3sin(^,.)))A j (21) 

hiP) = jj^, I (y^iP' co,{e,,))\ - /(X^(/5'sin(^,5)))A j (22) 

The first one, hs{P), is used to probe the response of the system to an external flux in the spatial 
planes (belonging to a 4d layer) while the second one, h^{l3), is used in a similar way for the 
planes containing the extra, transverse direction. 

(i) In the Strong phase (keeping (3 constant) the space-like helicity modulus vanishes (which 
is a clear signal of confinement); as we approach and eventually pass the phase boundary it 
becomes non-zero in the layer phase with a value that approaches 1 as /? increases further. On 
the other hand, the transverse h.m, hxiP ), remains zero throughout the transition since both 
phases exhibit confinement in the fifth direction. 

(ii) For the transition between the 5D Coulomb phase and the layer phase, hs{P) retains a 
value close to 1 for all values of P , since the four dimensional layers experience already a 
4d Coulomb-like phase, while hxiP ) vanishes for the layer phase; as the system crosses the 
critical point and enters the Coulomb phase it grows towards 1 as /3' increases further [|4l[T5l. 

3.2 The /?' = case 

On the axis defined by /3 = we consider the four-dimensional U(l) model. In this section our 
intention is to strengthen the arguments given in references [fT3| . We present numerical results 

11 



showing that we have a Coulomb phase only for T=0, in accordance with Fig. 12 of p3| . 
Our findings contradict the ones of ref. [! T2| that stipulates the existence of a Coulomb phase 
for Lt > 4 and P > Pc- The numerical results presented below show that we have a spatial 
confinement phase when the spatial lattice size Lg gets big enough, compared to the temporal 
size Lt {Ls >4:Lt). 

This behavior can be understood, following closely the ref. [fT3l . using simple theoretical 
arguments. In order to have hs{P) ~ 0, one must have at least two monopole loops (far apart) 
winding around the time direction with opposite orientations. A non contractible time-like 
monopole loop can, in principle, disorder all the spatial planes in the lattice. The probability to 
have one such loop passing through a given lattice site is, roughly, e~™'™"('^)^', where m„„„(/3) 
is the monopole mass. So the condition to achieve a probability of order 1 for a system to con- 
taining one (or two) wrapping monopole loops, is: 



L^ X e"'""™*^^^-^' ~ 1 



^ g^m„o„(/3)) 



L„ ~ e 3 



cf3 



(23) 



since m„„„(/?) is of order /?. Now, starting from the above equation we can make two state- 
ments. First, it predicts a pseudocritical coupling /3c ~ log(Ls) which was verified by our 
measurements and second, that as we go to smaller temperatures (bigger Lt), for sufficiently 
large Ls Eq. (23) is satisfied; hence the spatial planes becomes completely disorder and the 
Coulomb phase disappears. 

V= 2 X L.^ 



Ls=4 


L,=8 


L,=16 


Ls=24 



:;^U****^** 



P 



Figure 2: The mean value < \Pt\ > of the temporal Polyakov loop for /3 = and temporal size 
Lt = 2. < \Pt\ > goes to zero in the confining phase for Lg -^ oo. In the deconfinement phase 
(P > Pc{Lt)) < \Pt\ > is non-zero, approaching the value of one as /5 increases 

In Fig. 2 we show the mean value of the temporal Polyakov loop Pt{P) for Lt = 2. There 
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is an obvious continuous phase transition from the confining phase {jS < /3c — 0.90), where 
< I Pi I > is zero, to the deconfining phase where < \Pt\ > approaches the value of one. In the 
confining phase the free energy of a single static charge, relative to the vacuum, goes to infinity 
with Ls while it gets a positive value in the deconfining phase which vanishes as P increases. 
The mean value of the temporal Polyakov loop remains always non zero in the finite temperature 
phase as we switch on P' . This is the case presented in section 4. However this order parameter 
does not help us to characterize further the nature of the different phases. We arrived at the 
conclusion that the helicity modulus is a more promising quantity to study in detail the phase 
diagram. 





(a) 

Figure 3: The temporal helicity modulus ht for Lj = 2 (a) and Lt = 6 (b) and P = 0. Results 
from three different volumes are present. The value of ht increases with L^ in the deconfine 
region, for P bigger than a critical value, again with accordance with the scaling predictions of 
section (2.2.1). The transition for L^ = 2 is continuous, as opposed to the Lt = Q case, where 
we have a discontinuous behavior. 

In figure 3 we present our results for the temporal helicity modulus ht{P) for two different 
"temperatures" Lt = 2 and Lt = 6. The temporal h.m is zero in the confining phase for P < 
Pc{Lt) and non-zero for P > Pc{Lt) in the deconfining phase, indicating coulombic behavior. 
The signal for ht{P) in the deconfining phase (/3 > Pc) is being enhanced with increasing Lg, 
following the scaling relation ht ^ jf-. The transition point has only a weak dependence from 
the lattice volume showing convergence to a critical value Pc{Lt) with L^. We see that Pc{Lt) 
tends to smaller values as Lt decreases Q Another noticeable difference is the behavior of ht in 
the critical region. For Lt = 2 , ht{P) goes continuously to zero when P approaches Pc from 
above. For Lt = 6, on the other hand, the ht has an obvious discontinuity as P approaches Pc 



^(3c{Lt = 2) ^ 0.9008(3), /3c(it = 4) = 1.00340(1) and P^^Lt == 6) = 1.0094491(1). The results are taken 
fromref. d. 
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and the volume increases. This behavior indicates a different order for the phase transition, a 
second order phase transition for Lt = 2 and a first order for Lt = 6 (for details see ref. lfT3l ). 

In Fig. 4 we show the spatial helicity modulus hs{P) for Lt = 2 and spatial lattice sizes 
Lg = 4, 8, 16 and 24. The spatial helicity modulus is zero for /3 smaller than a critical value 
Pc{Ls) that depends strongly on Lg. We shall refer from now on to Pc{Ls) as the pseudocritical 
value, to distinguish it from the real critical value of P that comes from the temporal helicity 
modulus ht{P). For /? > Pc{Ls), hg takes non-zero values, increasing linearly as (3 takes bigger 
values . On the other side, the magnitude of this quantity decreases according to the ratio ~ ^, 
as we increase Lg and tends to zero for Lg — > oo. The pseudocritical value PdLg) increases 
with Lg as log(Ls) [13] and the ratio ^jjj\ tends to the value 0.5 for Lg > 16. In this way 
Pc{Lg) goes to infinity when Lg — i> oo. As a result the spatial h.m is always zero in the infinite 
volume limit and as a consequence we have a spatial confining phase. 
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Figure 4: The spatial helicity modulus hg for Lt = 2 and /3 = 0, versus the four dimensional 
coupling p. The pseudo critical value of P increases very fast, towards an infinite value, as 
the spatial lattice size Lg increases. The shift of the transition region to bigger values of P is 
obvious even in the smaller volumes. 



Finally we study the spatial helicity modulus for L^ = 4 and Lt = 6. The results are shown 
in Figs. 5(a)(Li = 4) and 5(h)(Lt = 6). We have, in general, the same situation as for Lt = 2. 
There is a pseudo-critical value of P that moves to larger values as Lg increases showing a 
strong dependence on Lt. The signal for Lt = 4 is clear only for L^ > 16. For Lt = 6 it seems 
that Lg = 16 is not enough but for Lg > 24 we get a clear displacement of Pc{Lg) to the right. 
In the region P < Pc{Lg) the spatial helicity modulus is zero (confining region). For P > Pc the 
spatial h.m scales with ~ ^ and tends to zero for L^ ^ oo. If we examine the Lt = 8 case for 
example, we would probably need volumes bigger than 8 x 32^ in order to get a clear picture 
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of the behavior of the system. From the previous observations we can say that hs{P) is zero for 
every value of P in the infinite volume limit, and consequently, we have spatial confinement for 
all temperatures different from zero. 

We conclude that the phase diagram on the /5, T = -^ plane has three phases: a confining 
phase for (3 < l3c{Lt), a temporal Coulomb - spatial confining phase for (3 > PdLt) U and the 
pure Coulomb phase for Lt ^ oo and (3 > jSc [fT3l[T4l . 
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Figure 5: In figures (a) and (b) we present the spatial h.m (hg) for /3 = and different temporal 
sizes, Lt = 4 and Lt = 6, for a variety of spatial volumes. The size of h^ decreases with 
Lg when P takes values bigger than the pseudo critical value, as predicted in section (2.2.1). 
The transition region moves clearly to the right as the volume increases, in agreement with the 
Lt = 2 behavior. 



3.3 Lt = 1 

In this case the temporal link is a Polyakov loop by itself and of course it is a gauge invariant 
quantity. The temporal plaquette becomes: 

0^,t{x) = d^{x) + 9t{x + ft) - 9^{x) - 9t{x) 

The two spatial links cancel each other, so in the U(l) case we get: 

Ofitix) = 6t{x + fi) - dt{x) 

The contribution to the action is: 



St = -j3 ^ cos(6'i(x + (j) - Ot{x)) 

z,l</x<3 



(24) 



^This phase is usually called deconfining phase. 
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The same applies for the "temporal-transverse" plaquettes and following the same steps as 
above we find that their contribution to the action is : 

S^' = -/?' J2 cosiOtix + 5) - Otix)) (25) 

From the above equations it can be observed that the temporal plaquettes decouple from the 
space and transverse ones. Equations (24) and (25) describe a 4D XY model with anisotropic 
couplings {P, P ). The three spatial links and the fifth, transverse link form a separate four di- 
mensional anisotropic U(l) gauge theory with two couplings, P and /5'. As a result the partition 
function of the model reduces to: 

■^(Lt=i) = ■^anisotropic 4D-xy ^ ■^anisotropic 4D-u(i) '^"^ 

and it describes two independent lattice field theories. 

The anisotropic AD — XY model, for /? =0, reduces to the three dimensional XY model 
which has a second order phase transition for /? = 0.4542 [1191 . The phase transition line 
continues to the (/?, /?') plane for smaller values of P as /?' increases and the critical value of P 
seems to tend asymptotically to the value of 0.1 as P' goes to infinity, u 

The 4D gauge model for /5 =0 reduces to a three dimensional U(l) gauge theory which 
is always in the confining phase. In the (/5, P ) plane we have a critical line which separates 
the strong confining phase from the four dimensional Coulomb phase. If we move along the 
diagonal, for example, where P = /5', we get the usual weak first order phase transition for 
p = p' = 1.001113 [|20l[Il. 

The above discussion can be summarized in the three dimensional plot of Fig. 6 . The 
vertical axis is for the temperature given in terms of the discrete variable Lt. The upper plane 
for Lt = 1 corresponds to "infinite" temperature while the lower plane for Lt = L^ = L^ 
corresponds to the zero temperature case. 

4 Study of the phase diagram for Lt = 2. 

In five dimensions the phase diagram at zero temperature is given in Fig. 1. [] For < P' < OAO 
and /? ~ 1. there is a critical horizontal line in the phase diagram separating the 5D strong 
confining phase from the layer phase. For P > 1 and P ~ 0.35 there is a critical vertical line 
that separates the layer from the 5D Coulomb phase. Our intention in this section is to explore 
the effects of finite temperature on our system and the most important, the feasibility (if any ) 
of a layer phase, through the study of the changes in the aforementioned phase line boundaries 



*For example, for the 4D XY model the critical value is at /? = /3 = 0.29(1) (see Fig. 6). 
'Lower plane (l/Lt ~ 0.0), see Fig. 6. 
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■•— 4d U(1) 
-•— 4dXY 
-♦— p =0 



Figure 6: Three dimensional phase diagram of the model, the vertical axis (l/Lt) represents 
the temperature. We present the critical curves for the three limiting cases. C4 and C5 are the 
four dimensional and the five dimensional Coulomb phases respectively. L stands for the layer 
phase at zero temperature. D^ is the temporal Coulomb - spatial confining phase for /3 = 0. 

and the phases themselves. To that end we move, first, on the line /3' = 0.20 in order to study 
the strong-layer phase transition at finite T; we know that for /5 =0 (subsection 3.2) there is 
phase transition for P ~ 0.90. Second, we move along the line P = 1.10, in order to study 
the layer-Coulomb phase transition at finite temperaturecl As we will explain in section 4.2 and 
using the Figs. 4 and 7 in order to have a clear picture of the behavior of the system for bigger 
values of P we need even bigger five dimensional volumes than those that we can presently 
achieve. 

Using the results presented in the two following sections we can argue that the layer phase 
disappears for Lt = 2 and becomes a deconfined phase with new properties which will be 
described below. We can also generalize the arguments and say that there is no layer phase in 
finite temperature for any temperature different from zero. The existence of the layer phase is 
based strongly on the existence of the Coulomb phase for /3 =0. However there is no Coulomb 
phase for /3' = at T 7^ as it is argued in ref. [fT3l . We also confirm this result (see subsection 
3.2). 



^We refer to the case of the plane {/3, (3 ) at l/Lt — 0.5 in Fig. 6. 
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4.1 Moving along the line /? =0.20 

We begin the investigation of the 5D anisotropic pure U(l) gauge model at finite temperature 
with what used to be called as a 5D strong-layer phase transition at zero temperature (Fig. 1). We 
utilize the helicity modulus hs{P), ht{P) in order to bring out the features of the transition and 
compare them with the T=0 and P = cases. As it is shown in Fig. 3 the first deviation from 
the zero temperature case comes from the fact that now, the transition line boundary between 
the two phases, is found at a lower value of /? = 0.90 in contrast with the value of /? = 1 .001 1 13 
for T = case. Another observation is that the values obtained here concerning ht{(3) are of 
the same order of magnitude as the ones, for the P = case; the only difference is the slight 
movement of the critical region to a value between 0.85 and 0.90. 

Moving now to a discussion of Fig. 7 and the spatial helicity modulus hs{P) we encounter 
many similarities with the results of subsection 3.2 : 

i) There is a pseudocritical value Pc{Ls) for each lattice size, with hg equal to zero for P < 
Pc{Ls), signal of spatial confinement. For P > Pc{Ls) the spatial helicity modulus hg increases 
with P, as one would expect from a Coulomb phase. But the transition point moves to higher 
and higher values of P as the spatial extent of the lattice {Lg) grows. What we see here is only 
a finite size effect that ceases to exist in the thermodynamic Lg ^ oo limit, 
ii) The magnitude of hg{P), calculated on a single 4d layer, decreases with Lg for the same 
value of P for P > Pc{Lg), following the ratio ~ j^. So we expect, as in the 4d case for 
P = 0, that the spatial helicity modulus tends to zero for all values of /? as L^ ^ oo (indicating 
spatial confinement); the phase transition to a Coulomb phase disappears together with the layer 
phase in the infinite volume limit. We mention also that the spatial-transverse helicity modulus 
(hgbiP)) remains zero throughout the transition. 
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Figure 7: The spatial helicity modulus hg is strictly zero for each volume V = 2 x L'^ until 
the pseudo-critical value Pc{Lg) is approached. For P > Pc{Lg) the hg tends to zero as -^ for 
constant p. 
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In Fig. 8(a) we present the temporal helicity modulus hf, also, in Fig. 8(b) we present 
the temporal-transverse helicity modulus ht^ versus P, for three different volumes. The two 
quantities have the same behavior: both take values equal to zero for P < 0.85 and non-zero 
for P > 0.85 and they increase with the lattice size Lg, indicating a coulombic behavior in the 
temporal direction. We also note that the Polyakov loop in the temporal direction, a result not 
shown here, is zero for P smaller than a critical value { Pc — 0.85) and tends to one for P > Pc. 
The transition, for the three quantities ht, ht5 and < \Pt\ >, concerning the strong confining 
phase (P < 0.85) to the deconfining phase (P > 0.85) is a continuous one . Although we do not 
analyze further the order of this phase transition we may guess that it may not be the case of a 
first order phase transition. 

All of the results obtained so far advocate to the disappearance of the layer phase at fi- 
nite temperature. The layer gives its place to a phase showing a confining behavior in the 4d 
subspaces (formed by the three spatial coordinates and the transverse one) and a coulombic 
behavior along the temporal direction. 



P =0.20 4d-layer 
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Figure 8: The temporal helicity modulus ht (a) and the temporal-transverse helicity modulus 
ht5 (b) for Lt = 2 and p' = 0.20 versus p. The ht is evaluated on the 4d-subspaces (Lt x L^) 
and scales as Lg for P > P^. The ht^ is evaluated on the whole lattice and scales as L^ for 

P>Pc. 



4.2 Moving along the line f] = 1.10 

As we have seen in the previous sections the system undergoes a continuous phase transition 
from the strong, confining phase, to a new phase. The transition point for /5 =0 is shown to be 
Pc ~ 0.90 and for P' = 0.20 it is slightly smaller being in the interval 0.85 < Pc < 0.90 region. 
In order to study the nature and the extent of the new phase, we choose to keep P fixed at the 
value of 1.10 and let P to vary. In Figs. 9(a) and 9(b) we present the spatial helicity modulus 
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hs{P ) and the spatial-transverse helicity modulus hs^i^P ) for three values of the volume. The 
hg and hg^ are zero, within the statistical error, for P' smaller than 0.445 signaling disordering 
in the spatial and transverse directions. This phase is the continuation of the P = phase to 
non zero values of P . The 3d U(l) theory obtained through dimensional reduction for P = 0, 
is extended (for < P < 0.445) to a 4d dimensionally reduced U(l) theory in the confining 
phase. We observe that the layer phase, consisting of a combination of 4d Coulomb phase and 
confinement in the extra dimension, becomes a deconfined phase. 

There is a critical region defined in the interval (0.445 < P < 0.450) in which a finite 
discontinuity in both quantities (hg, hg^) is shown up . For P > P^ the spatial helicity modulus 
is non zero and almost constant which is a characteristic of a Coulomb phase. The value of 
hsbiP') increases linearly with P' , following the lattice weak coupling expansion, approaching 
hs{P ) as P — > p. The values of hg and hsrj in Fig. 9 are divided by Lt and are independent of 
the spatial lattice size Lg. The spatial helicity modulus gives the renormalized coupling Pr = ^ 
of the 5D U(l) theory in the Coulomb phase which is fixed by the value of P = 1. 10 [4]. 

The temporal and the temporal-transverse helicity modulus ( not shown here), remain non 
zero and increase with P . Also the temporal Polyakov loop is non zero which is a signal of a 
finite temperature phase. 
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Figure 9: The spatial helicity modulus (a) and the spatial-transverse helicity modulus (b) as a 
function of P , measured for the temporal lattice size Lt = 2. The critical value of P remains 
constant with the lattice volume. 



By close inspection of Fig. 4 (/5 =0) and Fig. 7 (/? = 0.20), it becomes obvious that for 
a constant value of P the spatial helicity modulus is non-zero for some of the volumes that we 
used and it vanishes as the spatial volume increases beyond a definite value. For P = 1.10, for 
example, the lattice size Lg = 16 is enough to show the correct thermodynamic limit behavior. 
If we move to larger values of P, like P = 1.40, we have to use a spatial size of the order 
Lg > 24 in order to find the correct behavior. This is beyond our current computer capabilities. 
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In Fig. 10 we sketch, roughly, the phase diagram for L^ = 2 in the (P, P ) plane. There are 
three phases with different behavior of the observables we used: 

1. 5D confining phase with: Pt = 0, ht = 0, ht^ = hg^ = and hg = 

2. Finite temperature 5D Coulomb phase: Pt 7^ 0, ht,hs, ht5 and /i^s 7^ 

3. Dimensionally reduced 4d confining phase-temporal Coulomb Pt ^ 0, ht ^ ,ht5 7^ 
and hs, hs5 = 



CO. 
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0.94 

0. 
0.7 - 



P„h„h,5^ 



5D confining phase 

P„hs,h„hs5,h,5=0 



finite temperature 
5D Coulomb pliase 



Figure 10: A rough sketch for the phase diagram of the model for Lt = 2. There are three 
different phases, a 5D confining phase, a 5D Coulomb phase in finite temperature and a new 
one characterized as temporal Coulomb-4d confining. 

From the discussion in section 3.2 for /?' = we argue that the critical temperature for 
the appearance of the phase diagram of Fig. 10 it is the zero temperature. The reason is that 
the layer phase strongly depends on the existence of the phase transition in the Coulomb phase 
for P = 0. All the results we have presented in section 3.2 for T > and P = point to 
a 3d confining phase, in the infinite volume limit, for P larger than a critical value Pc{Lt)- A 
Coulomb phase does not seem to be the case. From this analysis we conclude that the phase 
diagram presented in Fig. 10 is reproduced for every temperature bigger than zero. Especially 
for P > Pc{Lt) and < /3 < PdLt) we have a 4d confining-temporal coulombic phase instead 
of a layer phase. Two charges are not anymore localized (confined) on a three dimensional 
subspace (brane) but the temperature gives the possibility of having interactions between the 
neighbor three dimensional subspaces. It seems that there are two characteristic correlation 
lengths in this deconfining phase. The correlation length given by the spatial string tension and 
a second one characterising the thickness of the brane given by the interaction in the transverse 
direction and the temperature. We did not study quantitatively these two correlation lengths at 
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finite temperature but we may easily see from the Fig. 6 what happens in two limiting cases. 
For Lt = I (infinite temperature) we have a 4d U(l) gauge theory in the strong confining phase 
and the two correlation lengths are indistinguishable and approach each other. For the zero 
temperature case on the other hand, there is no spatial string tension; we get a massless photon 
on the branes. Note that in this case the branes are characterized by zero thickness. In between 
these two limiting cases we expect a continuous change in the behavior depending strongly on 
the temperature. 

5 Discussion 

The extra dimensional models, like the brane models, are well studied mainly in the zero tem- 
perature case. But if we imagine that our brane world is a part of the Universe history then 
a study of the brane models at high temperature is required. In this paper we tried to do a 
first approach to this open problem, namely the behavior of brane models in high temperature 
(though neglecting the gravity effects). We believe that our toy model of five dimensional U(l) 
anisotropic lattice gauge theory has all the required essential characteristics. This model has a 
very rich phase diagram with respect to the temperature (see the discussion in Sections 3 and 4) 
summarized in Figures 6 and 10. 

Concluding we could note that the layer phase for zero temperature (with a massless photon 
on the brane and confinement in the extra dimensions) gives its place to a deconfined phase 
at non-zero temperature. In this phase the three spatial dimensions and the transverse one 
form a 4d subspace with confining properties, while the temporal direction shows a coulombic 
behavior. 
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